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Abstract 

We study the CP dependence of neutrino oscillation probability for all channels in arbitrary 
matter profile within three generations. We show that an oscillation probability for ~* i^fi 
can be written in the form P(i^e — i^fj.) = ^e/^ cos (J+Sg^ sin without any approximation 

using the CP phase S. This result holds not only in constant matter but also in arbitrary 
matter. Another probability for i/^ — > Vr can be written in the form P(f^ — > i^t) — cos6-\- 
B^r sin 6 + C^r + DfiT cos 26 + -B^t sin 26. The term which is proportional to sin 26 disappear, 
namely -B^r = 0, in symmetric matter. It means that the probability reduces to the same 
form as in constant matter. As for other channels, probabilities in arbitrary matter are at 
most the quadratic polynomials of sin (5 and cos 6 as in the above two channels. In symmetric 
matter, the oscillation probability for each channel reduces to the same form with respect to 
6 as that in constant matter. 



1 Introduction 

In solar and atmospheric neutrino experiments, the Ue deficit Q and the anomaly have 

been observed. These results strongly suggest the finite mixing angles 6*12 and 623 and the finite 

mass squared differences A 12 and A23, where Ajj = mf — rn^. Within the framework of three 

generations, there are two more parameters ^13 and 5 to be determined. About ^13, only upper 

bound is obtained from CHOOZ experiment Q and the information on the CP phase 5 is not 

obtained at all. In order to determine these parameters, several long baseline experiments using 

artificial neutrino beam will be planned and it is important to study the effect when the 

neutrino pass through the matter The main physics goal in these experiments is to measure 
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the value of 5. In this paper, we study the CP dependence of oscillation probability for all channels 
in arbitrary matter profile. 

Before giving our results, let us review the works on CP violation in three neutrino oscillation. 
At first we introduce the CP-odd asymmetry AP^^ = P{va ^p) — P{i^a — > i^/?)- In disappear- 
ance channel, AP^^ is exactly equal to in vacuum and independent of 5. However, AP^Jf is 
not always equal to in matter. This is due to the genuine CP violation and/or fake CP violation 
from matter effects. In the case of a = e, Kuo and Pantaleone Q have shown that P(fe i^e) 
does not depend on 6 in the context of solar neutrino problem. Therefore, AP^^ ^ arises from 
matter effects |^. However, in the case of a = ^, P(f^ v^^) has the CP-odd term in asymmetric 
matter profile as pointed out by Minakata and Watanabe 0]. We investigate the CP dependence 
in more detail in this paper. 

Let us consider the appearance channels. As the CP-odd asymmetry is proportional to sin 5 in 
vacuum, AP^^ ^ means that the discovery of CP violation. However, the situation completely 
changes when the matter effects are taken into account. Namely, AP^^ ^ does not always 
mean the existence of CP violation |^, because the fake CP violation due to matter effects exists. 
Here, it is difficult to separate genuine CP violation due to 5 from fake CP violation. One 
of the methods to solve these problems is to take into account mass hierarchy approximation 
IA21I ^ IA32I. Actually, some approximate formulae are given by Arafune et al. at low energy 
region and by Cervera et al. |jl^ and Freund |jll| at high energy region. 

Next, we introduce the T-odd asymmetry AP^ = P(z^q ^p) — P{^i3 ~^ t^o)- Krastev and 
Petcov [^] have shown that AP^ is proportional to sin (5 exactly in constant matter. Recently, 
Naumov |17] , Harrison and Scott [18| have derived the simple identity on the Jarlskog factor J 
|16| as A12A23A31J = A12A23A31J, where quantities with tilde represent those in matter. We 
can simply understand that AP^ is proportional to sin 5 from this identity. We have studied the 
matter enhancement of J |1£] taking advantage of this identity. Furthermore, Parke and Weiler 
have investigated the matter enhancement of the APg^ pO[| . 

There are some works on the deviation from constant matter. In long baseline experiments, we 
need to estimate the validity of constant density approximation because the earth matter density 
largely changes along to the path of neutrino. The matter profile of the earth is approximately 



expressed by Preliminary Reference Earth Model(PREM) Minakata and Nunokawa|22] 

give the oscillation probability using mass hierarchy and adiabatic approximations. For the 
distance less than L = 3000 km, the matter density fluctuation is small and the constant density 
approximation is valid. On the other hand, it has been shown that the fluctuation of the density 



cannot be ignored for the distance greater than L = 7000km |23, 24, 25, Furthermore, 
the constant density approximation is not valid in the case that the matter density profile is 
different from PREM and the asymmetric part exists. It is pointed out that AP^ has the term 
proportional to cos 5 in arbitrary matter. |27, 28|. We investigate this feature in more detail in 
this paper. 

In previous paper, we have proposed the new method applicable to constant matter. This 



^ Minakata and Watanabe have shown that Piue — » Ve) shghtly depends on 5 if we take into account the loop 
correction even in the standard model. In this paper, we do not consider the loop correction as these effects are 
safely neglected [M]. 
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method is to estimate the product of effective Maki-Nakagawa-Sakata matrix elements [[2 9| UaiU^i 
without directly calculating Uai- We have shown that the oscillation probability P{i'e — > t'^) is 
written in the linear combination |^ of cos 6 and sin 5 exactly |31] 



P(i/e i^fi) = cos 6 + Befj. sin 5 + Ce^, (1) 
in constant matter. In other channels, for example, 

P{^ti ^t) = Afj_r cos 6 + Bfj_r SUI 6 + C^r + DfiT COS 26. (2) 

It is found that the probability is quadratic polynomial of cos 6, sin 5, and the CP dependence 
was equal to in vacuum p2| . 

In this paper, we give the exact CP dependence of oscillation probability for all channels 
in arbitrary matter profile. For the purpose, we decompose the Hamiltonian H in the form 
H = {023Ts)H'{023ri)^ using 2-3 rotation matrix O23 and CP phase matrix = diag(l, 1,6**^). 
This decomposition plays a key role in our paper. As a result, we obtain the probability for 
f e f as 

P( z^e ~^ i^fj) — ^e/i cos 5 + B^,^ sin 5 + Ce^^. (3) 

This has the same form with respect to 5 as in eq.(|^) in constant matter. On the other hand, for 
t^fi ^ t^T, we show that the probability is given by 

P{j^fi i^t) = Afj_r cos 6 + Bfj^r siu 6 + Cf_ir + D COS 25 + E^r sin 25. (4) 

Comparing the equations (^) with (^), the probability in arbitrary matter profile has the term 
proportional to sin 25 which does not exist in the probability for constant matter. Furthermore, 
in the case of symmetric matter profile, we show that this additional term disappears, namely 
E^T = 0, and the probability reduces to the same form as in constant matter. 

2 CP Dependence in Arbitrary Matter Profile 

In this section, we study the exact CP dependence of neutrino oscillation probability in 
arbitrary matter profile. The Schrodinger equation for neutrino is 

i% = H., (5) 

where H is the Hamiltonian in matter and v is flavor eigenstate v = {Vf., v^, VrY' . We introduce 
the MNS matrix which relates the flavor eigenstate Va to the mass eigenstate Vi . The MNS matrix 
U in the standard parametrization is represented as 

U = 023TsOi3T\Oi2. (6) 



^ It is not so easy to obtain our result from the effective mixing and effective CP phase given by Zaglauer and 
Schwarzer iBOll 
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where Ts = diag(l, 1,6**^) and 



O 



23 



/ 1 \ 

C23 S23 
V -S23 C23 / 



(7) 



using the abbreviation Sij = sin 6ij and = cos 6ij. O13 and O12 represent 1-3 and 1-2 rotation 
matrix hke O23, respectively. By using this relation d^), we can rewrite the H as 



C/diag(0, A21, A3i)C/T + diag(a(t), 0, 0) 



2E 



023^6 Oi3Oi2diag(0,A2i,A3i)O?2O?3 + diag(a(t),0,0) TjO^a, 



(8) 



(9) 



where a{t) is matter potential defined by a{t) = 2\/^GFN{t)eE, and Gp, N{t)e, E are respectively 
Fermi constant, electron number density and neutrino energy. This equation (^) means that the 
Hamiltonian can be decomposed into two parts. One is 1-2 and 1-3 mixing part which contain 
matter effects. The other is 2-3 mixing and CP phase 8 part which does not contain matter 
effects. It is noted that this decomposition is garanteed by the relation 



O23r5diag(a(t),0,0)(O23r5)^ = diag(a(t),0,0). 

We can separate CP phase 5 from matter effects by taking advantage of this decomposition 
Changing v to v' as 

the Schrodinger equation d5[) is rewritten as 



(10) 



:ii) 



du' 



H'v' 



where 



2E 



Oi3Oi2diag(0, A12, Ai3)(Oi30i2)'^ + diag(a(t), 0, 0) 



(12) 



(13) 



We emphasise that the reduced Hamiltonian H' does not contain the 2-3 mixing and CP phase and 
is real symmetric. For anti-neutrino, we obtain the similar results by the replacements 6 — > —5 
and in a[t) — > —a{t). 

Next, we introduce the time evolution operator S{t) and S'{t) which is defined by the solution 
of the Schrodinger equation 

v{t) = S{t)v{^), v'{t) = S'{ty{^). (14) 
The relation between S{t) and S'{t) is determined by the transformation (111]) and is given by 



s{t) = (023r,)5'(t)(023r, 



(15) 
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By taking the component of (|T^, the relation between the time evolution operators for each 
flavour is given by 



See = S'ee, (16) 



S^e = S'C23 + S',,S2Ze'', (17) 



Sre = -S'f,eS23 + S'^^C23e'\ (18) 



S^^^, = S'^^C23 + S'^^C23S23e ' + S'^^C23S23e' + S'^^S23, (19) 
Sr>. = -S'^,C23S23 " S^^sige"*^ + S'^^cl^C^' + S',,C23S23, (20) 
Srr = S'A. - S'c23S23e-'' " 5' 023^236^' + S',A3^ (21) 



Here, Sap represents the transition amplitude for vp ^ Va- Se^, Ser and 5^,- are obtained from 
S/ie, Sre and Sri_i respectively by the replacements S'^^ — > 5^^, S —S. Substituting (p^-(^T|) 
into the relation 

P{l^a ^ 1^(3) = \Spa\^ (22) 

the oscillation probabilities in arbitrary matter profile are given by 



P{l^e- 


Ve) 


— Cee; 


(23) 


P{Ve- 




= Aefj, cos 5 + Befj. sin 5 + Ce^i, 


(24) 


P{Ve- 




= Aer COS 6 + Ber siu 6 + Cgr , 


(25) 






= A^^ COS 5 + Bfj,^ sin 6 + C^^ + cos 26 + E^^ sin 26, 


(26) 


P{Vr- 


Vr) 


= Arr COS 6 + B-TT sin 6 + Crr + Dtt COS 26 + E'tt sin 2(5, 


(27) 


P{^,^- 




= Afj^T COS 5 + B^T sin 5 + C^,- + D^^ cos 25 + E^j^t sin 25, 


(28) 



and the other probabilities Pii'^ t'e), -P(i^t ^e) and P{vr t'^) are obtained by the 
replacements S'^^ — > 5^^, and 6 ^ —6 in P{ve i^fi), Piye t^t) and P{i'^ Vt)-, respectively. 
Here all coefficients Ag^, ■ ■ ■ , E^j^r are constructed from the mixing angle 023 and 5' including 
matter effects. See appendix ^ for detail. The oscillation probabilities for "anti-neutrino" are 
also obtained by the replacements 6^—6 and a{t) — > —a{t). 

For these eqs.(p3|)-(p8|), we emphasize the following three points. First, the survival probability 
P(z/e — > i^e) in eq.(p3|) is compretely independent of CP phase 6. This coincides with the result 
by Kuo and Pantaleone |Q. Second, the transition probabilities P(fe Vfj) and P{i'e — > Vt) 
in eq.(^) and eq.(25) are linear polynomials of sin (5 and cos 6. These features coincide with 



the results in constant matter ||3l| , |32| . Third, ^(t'^ — > i^r), P{i^^i t^^i) and P{vr i^t) in 
eq.(p6[)-(p^) are at most quadratic polynomials of sin 6 and cos 6. 

We also comment the CP trajectory introduced by Minakata and Nunokawa. This is an orbit 



in the bi-probability space when 6 changes from to 27r |33, 34 1. Eq.(p^) shows that CP trajectory 
is exactly elliptic even in arbitrary matter profile. In addition, we point out that the dependence 
of 023 for the oscillation probabilities is completely understood from eqs.(p3D-(p7D. See appendix 
for detail. 

It is also noted that there are two features in asymmetric matter profile. First, the terms pro- 
portional to sin (5 and sin 2(5 are appeared in P^i'^ — > i'^) and P{i't — > ^t)- The term proportional 



5 



to sin 2(5 are appeared in Pii'/j, Vr) and P{vr v^). These terms do not exist in constant 
matter |31, 32]. Second, is not proportional to sin 5 in asymmetric matter as in constant 



matter [12|. In the next section, we describe these features in more detail. 



3 CP Dependence in Symmetric Matter Profile 

In this section, we study the CP dependence of P(z/q, up) in symmetric matter profile as 
special case of the previous section. In the case of symmetric matter along neutrino path, the 
time evolution operator 5" becomes symmetric matrix 

S'ali = S'pa^ (29) 

for flavour indices p8| , |35|| . As results, the relations between the coefficients of P{va T^e) 
eq.(p3[)-(p^) and P{ve z^o) are given by 

-^/xe — -^ein — B^^, C — C'g^, (30) 

See appendix P for detail calculation. The probabirity P{i'e ^a) have the same form with 



respect to 5 as eqs.(p3D-(25) in arbitrary matter proffie. 

On the other hand, applying the condition ( ^9|) to the probability (|2^ ) - (p8|) , we obtain the 
remarkable relations 

Bp,^ = Brr = E^^, = Err = E^r = E-r^ = 0, (32) 

where the detailed calculation is given in the appendix Using these relations, the oscillation 
probabilities (p6|)-(p^) have more simple form such as 

-P(t'/i I'fi) = cos 6 + C^^j + D^^ cos 25, (33) 

P{l>r Vt) = ^rr COS 5 + Crr + Drr COS 25, (34) 

P(f^ u^.) = Afj_r cos 5 + B^r sin 5 + C^^ + -D^r COS 25, (35) 

and P{ur i^fj.) is simply obtained by replacements Aa/3, B^p, Cap and Dap- Then, the coeffi- 
cients of P{i'r ^^i) are given by 

-^^T -^T^l B^f B^^, C^T C-ff^, DfMT Dt^i (36) 



where we use the condition ( ^91) for eqs.(|3q) or the unitarity for last equation. Here, the point is 
that the term proportinal to sin 25 is dropped in eq. ( p5| ) comparing with eq. (^) . The other point 
is that the terms proportional to sin 5 and sin 2(5 do not exist in P{v^ v^) and P{i'r ^t)-, 
As the results, the CP dependence of the oscillation probability for each channel in symmetric 



matter reduces to the same form as in constant matter |32|. This is the generalization of the 
result in our previous paper. 

Finally, we study the T-odd asymmetry AP^ = P(fa up) — P{'i^j3 — > I'a)- From the 
unitarity relation, we easily obtain 

APj^ = APj, = I^Pl, (37) 
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and in symmetric matter profile we obtain 

AP^^ = {Ae^-A^e)cos5 + {Bef,-Bf,,)sm6+{Ce^,-Cf,e) (38) 
= 2Bef,smS (39) 
= -4c23S23lm[5;;5;jsin5, (40) 

where we use the relations (pO|) in eg. (p9|) . In constant matter, Krastev and Petcov jl^ have 
pointed out that AP^^ is proportional to sin 6. Our result is applicable to the symmetric matter 
profile, which corresponds to the generalization of their result, even if the oscillation is non- 
adiabatic. 

Let us turn the case of arbitrary matter profile. in asymmetric matter is not proportional 

to sin J because the time evolution operator S' is not symmetric, namely 7^ S^^. More 
concretely speaking, the coefficients are not symmetric for flavor indices A^f^ 7^ Ap^, Bap 7^ —Bf^a, 
Cap 7^ Cpoi- We clarify the exact CP dependence of AP^^ in asymmetric matter although this 
fact is suggested using approximation |27, 

4 Summary 

We summarize the results obtained in this paper. We have studied the CP dependence of the 
oscillation probability P{i'a ^p) both in arbitrary and in symmetric matter profile. 

(i) In arbitrary matter profile, we have found that P{va ^p) is at most quadratic polinomial 
of sin 5 and cos 5. The CP dependences of the probabilities can be written as 

P{Ve ^ Ve) = Cee, (41) 
P{Va Vp) = COS 6 + Bap Sm 6 + Cap, 

for (a/3) = (en), (er), (^e), (re), (42) 
P{i'a — *■ i^p) = Aap COS 5 + Bap sin 5 + Cap + Dap COS 25 + Eap sin 25, 

for (a/3) = (/i/i), (fiT), (Tfi), (rr). (43) 

(ii) In symmetric matter profile, we have shown that the oscillation probabilities P(z^e — ^ i^x) 
have the same form as in arbitrary matter such as 

P(l/e ^ l^e) = Cee, (44) 

P{ve v^) = Ae^ COS 5 + B^^j. sin 5 + Ce^, (45) 

P(fe Vt) = Aer COS 5 + B^r siu 5 + Cgr • (46) 

Furthermore, we have shown that the CP dependences of other probabilities are written in 
the form as 

P{T^fi ^t) = Afj_r COS 5 + B^r sin 5 + C^r + D^r COS 25, (47) 

Pi^f, T^p) = COS 5 + Cf,f, + Df,^ COS 25, (48) 

P{Ur Vt) = COS 5 + Crr + Drr COS 25. (49) 

It is remarkable that the oscillation probability for each channel in symmetric matter reduces 
to the same form as in constant matter. 
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A Coefficients in Arbitraly Matter Profile 

In this appendix, we give exact CP and 2-3 mixing dependences of the oscillation probabilities 
in arbitrary matter profile. The probability for each channel is given by 

P{Ve^Ue) = Cee = Kf, (50) 
P{Ve^1^li) = Aefj,COs5 + Befj,S,Ul5 + Cepi, (51) 

A,^ = 2Re[5;;5;,]c23S23, (52) 



B,^ = -2lm[S';Xe]c23S23, (53) 

Ce/i = \S'f^e\ C23 + l^^el ^23' (54) 

P{Ve ^t) = ^er COS S + Bgr sin 6 + Cgr , (55) 

Aer = -2Re[S';X,]c23S23, (56) 

= 2Im[5;*,5;jc23S23, (57) 

Cer = I'S'liePsis + I'S'rePcIa, (58) 

P{i^ti^i^fi) = A^^cosS + Bfj_^s'mS + C^^ + Dfj_^cos26 + E^^sm26, (59) 

A,, = 2Re[(5;^4 + 5;,4)*(5;^ + 5;,)]c23S23, (eo) 

^M/x = -2Im[(5|,^c^3 + 5;^s^3)*(5';^-5|,^)]c23S23, (61) 

C/iAt = l'S'|i^pC23 + (I'S'^-rP + l'S'r^P)c23S23 + l'S'rrP*23 + "^^^[^ ixij.^'tt\'^23^23^ (*^2) 

D^, = 2Re[S':^S'^,]cl,sl,, (63) 

E^^ = 21ui[S';^S'^,]cl^sl„ (64) 

P{ur ^t) = Att- cos S + Brr sin 5 + C^r + Dtt COS 25 + Etj- sin 2(5, (65) 

A, , = -2Re[(5;^4 + 5;,ci3)*(S;^ + 5;j]c23S23, (66) 

= 2im[(5;^4 + s;,4)*(s;^-5;j]c23523, (6?) 

Crr = |S;^l'43 + (l'9;rl' + |5;Ml')44 + l'9;rl'4 + 2Re[S;;5;,]ci34, (68) 

Drr = 2Re[S'*^S'^,]c^,sl„ (69) 

Err = 2Im[5:;s;,]ci3S^3, (70) 

P(f ^ Ur) = A^r COS S + S^,- sin (5 + Cfj,r + -D/iT COS 26 + Efj^r sin 2(5, (71) 

V = -2Re[(5;^-5;j*(5;^ci3-5;,4)]c23S23, (72) 

B, r = 2Im[iS'^^-S'rrnS'r^4, + S'^rSl3)]c23S23, (73) 
C^tT = (I'S'^^P + I 'S'rr I ^)C23S23 + l'S'|trP*23 + l^Tixl'^ ^^23 ~ '^^^[^ i1ij,^'tt]'^23^23^ C''^) 
D,r = -2Rc[S'rlS'^r]43sl3, (75) 

E^r = -21m[S'r*^Slr]cl3sl„ (76) 

and the probabilities for the other channels P^i'n — ^ I'e), Piyr ^e) smd P{i'r ^ix) are obtained 

by the replacements S'^^ — S'j^^ and (5 — — (5 in P(fe — f;u), -P(fe — i^r) and -P(f/i — z^t), 



respectively. Prom these expressions, we can see that matter effects is renormalized in S'^p, which 
does not contain CP phase 5. Note that matter effects and CP effects are completely separated 
in the oscillation probability. The mixing angle ^23 is also separated from matter effects and all 
of the oscillation probabilities are quartet polynomials of C23 and S23. 

B Coefficients in Symmetric Matter Profile 

In this appendix, we give the relations of the coefficients of the probability in symmetric 
matter. We use the condition 

'S'a/3 = S'pa^ (77) 
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in symmetric matter profile. First, we calculate the relations between the coefficients of T- 
conjugate probabilities. From the oscillation probabilities (|5l| ) and (|7l| ) in appendix |^ and the 
symmetry of S' (|77|), we obtain 



A,^ - A^, = Ac23S23Re[S';Xe - S'elS',,] = 0, (78) 

Be^ + B^e = -4c23S23lm[5;:5;, - = 0, (79) 

Cefi — C^e = (|5'^el^ ~ I 'S'e^i I ^ ) ^23 + (I'^TePl — l'S'eTp)'523 = 0) (^0) 

V - Ar^. = 4Re[(S;^ - S;j*(S;^ - S'^^)] = 0, (81) 

B^r + B^^ = -4Im[(5;^ - S'^^nS',^ - S'^,)] = 0, (82) 

C,,r - = {\S',,\' - K/)isi, + 4) = 0, (83) 

D^r - D^^ = -4(|5;,|2 - \S',/)4,sl, = 0. (84) 

The relations between the coefficients for I'r are obtained in the same way. 

Second, we calculate the coefficients of sin 5 and sin 2(5 in Pii'^ — > f^) and P(z^r t^t) and the 
coefficients of sin 25 in P{i'^ v^) and P{vr — > t^t)- From the concrete expression of oscillation 
probabilities (p9[), ( |65[ ) and ( |7l[ ) in appendix^, we obtain 

B^^ = -2c23S23lm[(5;^4 + S',Az)*{S',>. - S'^r)] = 0, (85) 

E^^ = 2lm[S':^S'^M,sl, = 2lm[\Sl,,\^]cl,sl, = 0, (86) 

Brr = 2c23S23lin[{S'^^sl, + S',,4,r{S',^ - S'^,)] = 0, (87) 

Err = 2lm[S';^S'^r]4,sl, = 2lm[\S'^r\']4,sl, = 0, (88) 

E^r = -2lm[S';^S'^r]4,sl, = -2lm[\S'^r\^]4,sl, = 0, (89) 

Er^ = 21m[S';rS'r^]4,sl, = 2lm[\S'r/]43sh = 0, (90) 



from the condition (^). 
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